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ABSTRACT 

The broad band spectra of two Swift/BAT AGNs obtained from Suzaku follow-up observations are 
studied: NGC 612 and NGC 3081. Fitting with standard models, we find that both sources show 
similar spectra characterized by a heavy absorption with Ah — 10 24 cm -2 , the fraction of scattered 
light is / sca t = 0.5 — 0.8%, and the solid angle of the reflection component is VL/2tt = 0.4 — 1.1. 
To investigate the geometry of th e torus, we apply numerical spectral models utilizing Monte Carlo 
simulations by llkeda"et al. (2009) to the Suzaku spectra. We find our data are well explained by this 
torus model, which has four geometrical parameters. The fit results suggest that NGC 612 has the 
torus half opening-angle of ~ 60° — 70° and is observed from a nearly edge-on angle with a small 
amount of scattering gas, while NGC 3081 has a very small opening angle ~ 1 5° and is observe d 
on a face-on geometry, more like the deeply buried "new type" AGNs found bv lUeda et all (pOOl . 
We demonstrate the potential power of direct application of such numerical simulations to the high 
quality broad band spectra to unveil the inner structure of AGNs. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

The strong correlation between the mass of a super- 
mass ive black hole (SMBH) and that of the galactic bulge 
(e.g.. lMagorrian et al.lfl998l iMarconi fc Huntll2003f) sug- 
gests a fundamental link between the growth of SMBH 
and galaxy evolution. Theoretical models predict that 
most SMBHs in galaxies experien ce a heavily obscure d 
phase in their growth stage (e.g., iHopkins et al J [2005) . 
Indeed, studies based on population synthesis models of 
the cosmic X-ray background (CXB) suggest that heav- 
ily obscured AGNs, whose line-of-sight hydrogen column 
density (Ah) is greater than 10 23 5 c m" 2 , are a signifi- 
cant fraction of the AGN population (|Ueda et al.l 120031 : 
iGilli et alj|2007ft . Due to the difficulty of detecting them 
in most energy bands, however, our understanding of 
heavily obscured AGNs (including "Compton-thick" ones 
with Ah > 10 24 cm~ 2 ) is very scarce even in the local 
universe. 

Sensitive hard X-ray observations above 10 keV, where 
the penetrating power overwhelms photo-electric ab- 
sorption, provide fruitful information about this pop- 
ulation, except for heavily Compton thick (Ah > 
lO 24 5 cm" 2 ) objects. Recent all sky hard X-ray surveys 
perfo rmed with Swift/BAT (15-20 keV: iTueller et al.l 
20081) and INTEGRA L (10-100 keV: lBassani et al.ll2006l: 
Krivonos et al.l |2007|) are ideal for this purpose with 
much less selection biases than surveys at lower energies. 

Our team have been working on a systematic follow-up 
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observation program with Suzaku of Swift/BAT detected 
AGNs whose broad band X-ray spectra were poorly (or 
never) studi ed previous l y, targ eting obscured objects in 
most cases. lUeda et ail (|2007t ) discovered deeply buried 
AGNs that exhibit very small fractions of scattered soft 
X-rays (< 0.5%) with respect to the transmitted compo- 
nent, with strong reflection signals most probably coming 
from the inner wall of the Compton-thick tori. Further 
studies of six Swift AGNs bv lEeuchi etail ((20091) (Paper 
I hereafter) show that they could be classified into two 
types, "new type" AGNs with a small scattering frac- 
tion and strong reflection strength, and "classical type" 
ones with a larger scattering fraction and weaker reflec- 
tion. These types are consistent with SMBHs surrounded 
by geometrically thick and thin t ori, respectively . Using 
an INTEGRAL selected sample, iComastri et~atl (|2009h 
also suggest that there are distinct AGN populations 
of new and classical types, although the result depends 
on whether the absorption for the reflection compo- 
nent is considered o r not in the spectral model (see 
IComastri et al.ll2010l ). Due to the limited number of ob- 
jects in the sample studied we are far from reaching a 
consensus on the torus structure and its dependence on 
various parameters like the AGN luminosity, Eddington 
ratio, and properties of the host-galaxy, for the whole 
AGN populations. 

High quality broad band X-ray spectra give unique in- 
sight into the structure and geometry of the central re- 
gion of AGNs. Most previous studies, however, relied 
on phenomenological spectral models where the detailed 
geometry of the torus is not taken into account; usu- 
ally, an analytical formula for the Compton reflection 
from matter with infinite optical depths is simply as- 
sumed for the reprocessed emission, and absorption col- 
umn density of the transmitted component is treated in- 
dependently. Monte Carlo simulation is a powerful tool 
to reproduce realistic spectra from AGNs with a com- 
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Table 1 

List of Targets 



SWIFT Optical/IR Identification R.A. (J2000) Dec. (J2000) Redshift Classification 



J0I34.I-3625 NGC 612 01 33 57.74 -36 29 35.7 0.0298 Seyfert 2 

J0959.5-2258 NGC 3081 09 59 29.54 -22 49 34.6 0.0080 Seyfert 2 



Note. — The position, redshift, and classification for each source is taken from the NASA/IPAC Extra- 
galactic Database. 



plex structure of the torus, which may not always have 
a sufficiently large optica l depth for Compton scatter- 
ing. Recently, llkeda et "all (|2009f ) have developed such a 
Monte Carlo code that can be applicable to the broad 
band X-ray spectra with several free parameters describ- 
ing th e torus geometry. Similarly, iMurphv fc Yaaoobl 
(2009) also studied the numerical spectra from a toroidal 
torus, known as the MYTORUS mode0, although we do 
not adopt this model here because the opening angle of 
the torus is fixed at 60°. Applying such models directly 
to the observed spectra, we can obtain more accurate 
constraints on the inner structure of AGNs than from 
the standard previous analysis. 

In this paper, we present the results of detailed X- 
ray spectral analysis of Suzaku data of two Swi/f/BAT 
AGNs, Swift J0134.1-3625 (NGC 612; z = 0.0298) and 
Swift J0959. 5-2258 (NGC 3081; z = 0.0080), whose si- 
multaneous broad band spectra were not available before. 
NGC 612 is a powerful radio galaxy, which was origi- 
nally classified as Fanaroff-Riley ( FR; |Fa naroff fc Rilevl 
11974 II type by iMorganti etafl (|1993t l This object 
hosts prominent double radio sources; the eastern lobe 
has a bright hot spot near its outer edge, while the 
western one has a jet-like structure. Since the for- 
mer and latter morphology cor respond to those of the 
FR I and II types, respectively, iGopal-Krishna fc Wiital 
( 2000) classifies it as a hybrid morphology radio source. 
NGC 612 shows an optical spectrum of Seyfert 2 galax- 
ies but the intensity of the [O III] emission is very weak 
(jParisi et all 12009ft . IWinter et all ()2008[ ) present the X- 
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ray spectrum observed with XMM-Newton, obtaining a 



large hydrogen column density of Nn 
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with an apparently very flat power-law index of ~0.3, 
suggestive of a reflection-dominant spectrum below 10 
keV. NGC 3081 is a Seyfert 2 hosted by a barred 
galaxy. This ob ject has three rings associated with a 
tidal interaction ( Freeman et al.ll2000|) . and e xhibits time 
varia ble polarization in the optical band (jJoshi et al.l 
1989). Th e strong [O III] emission is observed from 
this ob ject (jStorchi-Bergmann et al.lH995[ ). iMoran et al.l 
(2001) report a large absorption column density (JVh — 
10 23 ' 7 cm~ 2 ) from the X-ray spectrum observed with 
ASCA. 

Section [2] describes the observations and data reduc- 
tion of the two sources. We first analyze the spectra 
with standard spectral models in Section [3J and then 
present the results of application of the torus model uti- 
lizing Monte Carlo calculation by llkeda et al.l (|2009f) in 
Section 2J The implications of our results arc discussed 
in Section [5] We adopt the cosmological parameters (H , 
n m , Q x ) = (70 km s _1 Mpc" 1 , 0.3, 0.7) throughout the 



2. OBSERVATION AND DATA REDUCTION 
2.1. Observation 



We observed NGC 612 and NGC 3081 with Suzaku in 
2008 May and June, respectively. The basic informa- 
tion for our targets is summarized in Table [TJ Suzaku 
(jMitsuda et al.l 120071 ) carries four X-ray CCD cameras 
called the X-ray Imaging Spectrometer (XIS-0, XIS-1, 
XIS-2, and XIS-3) as focal plane imager of four X-ray 
telescopes, and a non-imaging instrument called the Hard 
X-ray Detector (HXD) consisting of Si PIN photo-diodes 
and GSO scintillation counters. XIS-0, XIS-2, and XIS-3 
are front-side illuminated CCDs (FI-XISs), while XIS-1 
is the back-side illuminated one (BI-XIS). To maximize 
the effective area of the HXD, the targets were observed 
at the HXD nominal position, which is about 5 arcmin 
off-axis from the averaged optical axis of the XISs0 

We analyze only the data of the XISs and the 
HXD/PIN, which covers the energy band of 0.2-12 keV 
and 10-60 keV, respectively The fluxes above 50 keV are 
too faint to be detected with HXD/GSO. Table [U shows 
the log of the observations. The net exposure of each 
target is about 45 ks. Bec ause XIS-2 became unoperat- 
able on 2007 November 7 (jDotani et al.ll2007t ). no XIS-2 
data are available for both objects. For the XIS obser- 
vations, we applied spaced-ro w charge injection (SC I) to 
improve the energy resolution (jNakaiima et al.|[2008|) ; for 
instance, it reduces the full width at half maximum of the 
55 Fe c alibration source fr om ~ 230 eV to ~ 160 eV for 
XIS-0 (jOzawa et al. 200j|). To constrain the broad band 
spectra above 60 keV, we also utilize the Swift/BAT spec- 
tra covering the 15-200 keV band, integrated over the 
first 22-months of Swift operations. 

2.2. Data Reduction 

The Suzaku data are analyzed by using HEAsoft ver- 
sion 6.7 and the latest version of CALDB on 2009 De- 
cember 3. For the XIS data, we analyze the version 2.2 
cleaned events distributed by the Suzaku pipeline pro- 
cessing team. In extraction of the light curves and spec- 
tra, we set the source region as a circle around the de- 
tected position with a radius of 1.5 arcmin, where about 
75% of the total source photons are accumulated, to max- 
imize the signal-to-noise ratio. The background for the 
XIS data is taken from a source-free region in the field of 
view with an approximately same offset angle from the 
optical axis as the source. For the non X-ray background 
of the HXD/PIN data, we use the so-called "tuned" back- 
ground model provided by the HXD team. Its system- 
atic errors are estimated to be ~ 0.97% at a lcr confi- 
dence level in the 15-40 keV band for a 40 ks exposure 
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Table 2 

Observation Log 



Target 


Start Time (UT) 


End Time 


Exposure a (XIS) 


Exposure (HXD/PIN) 


SCI b 


NGC 612 


2008 May 20 16:19 


May 21 20:08 


48.5 ks 


41.3 ks 


On 


NGC 3081 


2008 Jun 18 21:49 


Jun 19 19:33 


43.7 ks 


42.4 ks 


On 



a Based on the good time interval for XIS-0. 

With/without the spaced-row charge injection for the XIS (Nakaiima et al.l l20Q8). 



nlime: 5760 s 



(a-1) NGC 612 XIS 



+ 



reduced f (dof) - 0.62 (17) 



(a-2) NGC 612 HXD/PIN 



Bin time: 5760 



reduced t (d af) = 1.33 (17) 



2x10* 4x10 4 6x10 4 8x10 4 10 E 
Time (s) 

Start Time 14606 17:08:03:795 Stop Time 14607 20:20:03:795 



2x10 4 4x10 4 6x10 4 8x10 4 1 
Time (s) 

Start Time 14606 17:09:56:169 Stop Time 14607 20:21:56:1 



Bin time: 5760 



Bin time: 5760 



(b-1) NGC 3081 XIS 



++ 



++ 



+++++ + 

+ ++ + 



reduced x 2 (dof) = 5.22 (13) 



(b-2) NGC 3081 HXD/PIN 



t 



reduced x 2 (dof) = 1.81 (13) 



2x10* 4*10« 6x10 4 8x10" 10 B 
Time (s) 

Start Time 14635 23:05:33:892 Stop Time 14636 19:53:33:892 



2x10' 4x10" 6x10' 8x10" 10 B 
Time (s) 

Start Time 14635 22:58:33:184 Stop Time 14636 19:46:33:184 



Figure 1. The background subtracted light curves of Suzaku. One bin corresponds to 96 minutes. The numbers listed in each panel 
represent the value of reduced x 2 with the degrees of freedom for the constant flux hypothesis. Left: The light curves of the XIS in the 
2-10 keV band. The data from the XIS-0 and XIS-3 are summed. Right: The light curves of the HXD/PIN in the 15-40 keV band. 



(jMizuno et al.| [2008). Since our exposures are ~ 40 ks or 
longer, we expect that the error is even smaller than this 
value. The CXB spectrum simulated with the HXD/PIN 
response for a uniformly extended emission is added to 
the non X-ray background spectrum. 

2.3. Light Curves 

Figure [TJ shows the background-subtracted light curves 
of our targets obtained with the XIS and HXD/PIN in 
the 2-10 keV and 15-40 keV bands, respectively. To min- 
imize any systematic uncertainties caused by the orbital 
change of satellite, we merge data taken during one or- 
bit (~ 96 minutes) into one bin. Then, to check if there 
are any significant time variability during the observa- 
tions, we perform a simple y 2 test to each light curve 
assuming a null hypothesis of a constant flux. The re- 
sultant reduced \ 2 value and the degrees of freedom are 
shown in each panel. As noticed from Figure [IJ the 2-10 
keV flux of NGC 3081 increased by a factor of 1.5 af- 



ter ~ 20 ks from the start of the observation. Then, a 
flux decline is suggested between ~ 50 ks and ~ 60 ks 
particularly in the 15-40 keV band. Thus, we divide the 
observation of NGC 3081 into three different time re- 
gions, 0-20 ks (Epoch 1), 20-60 ks (Epoch 2), and 60-80 
ks (Epoch 3) measured from the observation start. By 
contrast, no significant time variability on a time scale 
of hours are detected from NGC 612. Hence, we analyze 
the time-averaged spectra over the whole observation for 
NGC 612. 

2.4. BAT Spectra 

It is known that the incident photon spectra of Seyfert 
galaxies are roughly approximated by a power law with 
an exponential cutoff (cutoff power- law model), repre- 
sented as AE~ T exp (— E/E cut ), where A, T, E cut are 
the normalization at 1 keV, photon index, and cutoff 
energy, respectively. We analyze the Swift/BAT spec- 
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(a) NGC612 



Table 3 

Cutoff Energies (_E C ut) determined 
by the BAT spectra 



Q/2tt 


NGC 612 


NGC 3081 





> 315 


> 338 


x 7d.o.f. 


16.4/6 


11.6/6 


2 


> 319 


> 293 


X 2 /d-o.f. 


10.3/6 


5.8/6 



Note. — The unit of E c 



koV. 



tra in the 15-200 keV band to constrain E cut . Here we 
take into account possible contribution from a Compton 
reflection component fro m optically thick, cold matter , 
utilizing the pexrav code lMagdziarz fe Zdziarskil {1995). 
The relative intensity of the reflection component to that 
of the intrinsic cutoff power-law component is defined as 
R = fl/2n, where f2 is the solid angle of the reflector 
(R = 1 corresponds to the reflection from a semi-infinite 
plane) . 

In the analysis of the Swift/BAT spectra, we assume 
R = or 2 as the two extreme cases just to evaluate the 
effects of including the reflection components, as done in 
Paper I. The inclination angle is fixed at 60°. To avoid 
strong coupling between the power-law slope and cutoff 
energy, we fix th e photon index at 1 . 9, the canonical slope 
for AGNs (e.g., iNandra fc Poundslfl99l . Table H gives 
the fitting results for E cnt ', we find that E cnt is greater 
than ~ 300 keV for both targetsQ Accordingly, we fix 
it at 300 keV (or 360 keV for consistency with the Ikeda 
model) in the following spectral analysis. 

3. ANALYTICAL MODELS 

We perform the spectral fitting to the Suzaku data in 
the same manner as Paper I. We start with the sim- 
plest model for each target, and if we find that the 
fit with a simple model does not give a physically self- 
consistent picture or that the fit is significantly improved 
by introducing additional parameters, then we adopt 
more complicated models. We use only the Suzaku XIS 
and HXD/PIN data throughout this stage, and finally 
perform the simultaneous fit of XIS, HXD/PIN, and 
Swift/BAT spectra with the selected model to obtain the 
best-fit parameters. 

The spectra of FI-XISs are summed, and the rela- 
tive normalization between the FI-XISs and the PIN 
is fi xed at 1.18 based o n the calibration of Crab Neb- 
ula (|Maeda et al.1 12008T) . Those of BI-XIS and BAT 
against FI-XISs are set as free parameters. The Galac- 
tic absorption (iV^ al ) is always included in the mod- 
els, whose hydrogen column d ensity is fixed at va l- 
ues obtained from the H I map Kal berla et al.l ((2005), 
available with the nh program in the HE As oft pack- 
ag e. We adopt the photoelectric absor p tion c ross section 
bv iBalucinska-Church fc McCammonl (Jl992) ("b cmc ). 
Different from Paper I, we allow the iron abundance to 
be a free parameter by us ing zvphabs because non So- 
lar values (as defined by lAnders fc Grevessel 119891 ) are 
required to explain the Suzaku spectra, while Solar abun- 

7 This conclusion is unchanged when we fix the inclination angle 
at 30° or 80° , indicating that E cu t is not sensitive to the assumed 
inclination angle. 
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Figure 2. The confidence map in A% 2 (color scale) with respect 
to the strength of the reflection component (R = Q/2n) and its ab- 
sorption (N™ a ) for NGC 612 (top) and NGC 3081 (bottom). The 
dashed and solid curves correspond to the ltr and 2a confidence 
level for two interesting parameters, respectively. 



dances are adopted for the other metals throughout our 
analysis. 

We use the same three models as defined in Paper I, 
consisting of an absorbed transmitted component, a scat- 
tered component, and/or an absorbed reflection compo- 
nent, with an iron-K emission line: 

• Model A: transmission + scattering + iron line0 

• Model B: transmission + scattering + iron line + 
absorbed reflection^ 



• Model C: transmission with dual absorber 
tering + iron line + absorbed reflection^ 



seat- 



In our analysis, we adopt an unabsorbed power law with 
the same photon index as the incident continuum to de- 
scribe the scattered component, ignoring any emission 
lines from the photo-ionized gas. Note that here we only 
introduce a single absorber for the reflection component 
as the simplest approximation, although we expect both 
absorbed and unabsorbed ones from the torus as well as 
that from the accretion disk, as described in next Sec- 
tion. The pexrav component in each spectral model 
(see footnote) represents only the reflection component 
not including the direct one by setting R < 0, and the 

8 In XSPEC nomenclature, zvphabs*zhighect*zpowerlw + 
const*zhighect*zpowerlw -j- zgauss 

9 In XSPEC nomenclature, vzphabs*zhighect*zpowerlw + 
const*zhighect*powerlw 4- zgauss + zvphabs*pexrav 

10 In XSPEC nomenclature, zvphabs*zpcfabs*zhighect*zpowerlw 
+ const*zhighect*zpowerlw + zgauss + zvphabs*pexrav 
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inclination angle is fixed at 60°. Theoretically, the equiv- 
alent width (E.W.) of the iron-K emission line with re- 
spect to the re flection componen t, E.W. rcfl , is expected 
to be ~1 keV (jMatt et al.l ll991j) . Since this value de- 
pends on the geometry of the reflector as well as the iron 
abundance, we regard the result as physically valid if 
E.W. rcfl = 0.5-2 keV. No other emission lines than iron 
Ka are significantly detected from the spectra. Consid- 
ering calibration uncertainties in the energy response of 
the XISs, we fix the la line width of the iron-Ka emis- 
sion at the averaged value of the (apparent) line width 
of the 55 Fe calibration source at 5.9 keV: 45 eV and 47 
eV for NGC 612 and NGC 3081, respectively. 

3.1. NGC 612 

Model B is adopted as the most appropriate model of 
NGC 612. We obtain {\ 2 ,v) = (90.8,85) with Model A 
and (x 2 ,^) = (86.9,84) with Model B from the Suzaku 
spectra, where v is the degree of freedom. Thus, the im- 
provement of the fit by adding a reflection component 
is found to be significant at 94% confidence level by an 
F-tcst. No significant improvement is found with Model 
C. For this target, the absorption to the reflection com- 
ponent A H ofl is linked to that for the transmitted compo- 
nent, because making them independent does not give a 
better fit over the statistics (see below). Since the E.W. 
of the iron-K line with respect to the reflection compo- 
nent is E.W. refl = 0.7 ± 0.1 keV, the model is physically 
self-consistent; the iron abundance of NGC 612 is roughly 
half of the Solar value, 0.54±g;^, obtained from the si- 
multaneous fit of the Suzaku and BAT spectra. 

To examine the degeneracy in the fitting parameters, 
m Figure [2] (top) we show the confidence map (in terms 
of Ax 2 ) with respect to the strength of the reflection 
component (R) and its absorption (Ay cfl ), based on the 
Model B fit including the BAT data. Here we do not link 
A H to that of the transmitted component (Ah), and 
explore a region of A|j ofl < 1.2 x 10 24 cm" 2 , the upper 
limit obtained for Ah- The contours give the condidence 
levels at la and 2a for two interesting parameters. As 
noticed, while a wide range of A H (> 6 x 10 23 cm.- 2 , la) 
is allowed, we can constrain the reflection strength to be 
R ~ 0.6 for A£ fl < 1.1 x 10 24 cm" 2 and R ~ 0.3 - 0.7 
otherwise. The case of an unabsorbed reflection compo- 
nent (A H efl = 0) or no reflection component (R = 0) is 
rejected at > 99% confidence level, which corresponds to 
Ax 2 = 9.21. 

3.2. NGC 3081 

First, we analyze the Suzaku spectra integrated over 
Epoch 1. We obtain (x 2 ,^) = (120.1,61) with Model 
A and (x 2 ,^) = (80.2,59) with Model B, and thus the 
improvement of x 2 is significant at > 99% confidence 
level by an F-test. No significant improvement is found 
with Model C. Finally, since positive residuals remain in 
the energy band below 1 keV, we add the vape<0 in 
XSPEC, a spectral model from an optically-thin thermal 
plasma, whose iron abundance is linked to that in the 
absorber of the transmitted component. This yields a 
further significantly better fit at > 99% confidence level 

11 http: //exc .harvard. edu/atomdb/ 



with (x 2 ,z/) = (50.7,57). The E.W. of the iron-K line 
with respect to the reflection component is E.W. lefl = 
1.5 ± 0.3 keV, which is self-consistent. Thus, we adopt 
Model B+vapec as the best fit model for Epoch 1. 

By fitting the Epoch 2 and 3 spectra of Suzaku 
with the same model, we obtain acceptable fits with 
( x 2 ,v) = (139.7,114) and (19.1,30), respectively. The 
best-fit parameters in the three epochs are plotted in 
Figure [3] except for the iron abundance and those of the 
thin thermal component, which are expected to show 
no time variability. We find only a weak indication 
that the column density changed during the observa- 
tion from A H = (113 ± 10) x 10 22 cmT 2 (Epoch 1), 
to A H = (98 ± 6) x 10 22 cur 2 (Epoch 2), and then 
A H = (86 ± 10) x 10 22 cm~ 2 (Epoch 3) in addition 
to the unabsorbed power-law luminosity, which varied 
from £ 2 _i = (3.0 ± 0.2) x 10 42 ergss" 1 (Epoch 1), to 
L 2 _io = 3.8±o 2 x 10 42 ergss" 1 (Epoch 2), and L 2 -io = 
2 - 7 ±8J x 10 42 ergss" 1 (Epoch 3) in the 2-10 keV band. 
The other parameters are found to be consistent with be- 
ing constant among the three epochs within the errors. 
The significance of the variability of the column density 
is marginal, as the null hypothesis probability of a con- 
stant value is found to be 15% from a x 2 test. Thus, we 
sum the Epoch 1, 2 and 3 data of Suzaku and discuss the 
time averaged spectra in the following analysis. 

To best constrain the spectral parameters of 
NGC 3081, we perform the simultaneous fit to the 
time-averaged Suzaku and BAT spectra with the Model 
B+vapec model. This yields (xV) = 0(199.28,199) 
and E.W. rcfl = 1.0 ± 0.1 keV, and thus is physically self- 
consistent. The iron abundance with respect to Solar is 
0.89 ± 0.07, and the temperature of the plasma is found 
to be kT = 0.26 ± 0.02 keV with an emission measure of 
n 2 V ~ 1.5 x 10 63 cm- 3 . 

Figure [2] (bottom) shows the confidence contour map 
with respect to R and based on Model B (includ- 

ing the BAT data) for NGC 3081. Unlike the case of 
NGC 612, the solution is well constrained (R ~ 0.8 — 1.0, 
la) and we do not sec strong degeneracy in the fitting 
parameters. Again, neither the case of an unabsorbed 
reflection component (A|j cfl — 0), nor no reflection com- 
ponent (R — 0) is allowed at > 99% confidence level. 

3.3. Results Summary of Analytical Models 

We summarize the best-fit models and parameters in 
Table |U The observed fluxes in the 2-10 keV and 10-50 
keV bands, and the estimated 2-10 keV intrinsic lumi- 
nosities corrected for absorption are also listed. Figure [4] 
(left) shows the observed spectra of the FI-XIS (black), 
the BI-XIS (red), and the HXD/PIN (magenta) folded 
with the detector response in units of counts s _1 keV -1 , 
together with the unfolded BAT spectra (blue) in units 
of photons cm- 2 ks -1 keV -1 . The best-fit models are 
superposed by solid lines. In the lower panels, the cor- 
responding data-to-model residuals in units of x (i- e -j 
normalized by the la statistical error in each bin) are 
plotted. Figured] (right) shows the best- fit spectral mod- 
els in units of EFe without Galactic absorption, where 
the contribution of each component is plotted separately; 
the black, red, blue, cyan, magenta curves correspond to 
the total, transmitted component, reflection component, 
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Figure 3. Time variability of the best-fit parameters of NGC 3081 obtained with the analytical model: from left to right and top to 
bottom, (a) the line-of-sight hydrogen column density for the transmitted component, (b) the power-law photon index, (c) the fraction of 
the scattered component relative to the intrinsic power law, (d) the center energy of the iron-K emission line at the rest frame, (e) the 
equivalent width of the iron-K line with respect to the whole continuum, (f) the line-of-sight hydrogen column density for the reflection 
component, (g) the solid angle of the reflection component, (h) the 2-10 keV intrinsic luminosity corrected for the absorption. 
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Figure 4. The observed spectra (left) and the best-fit spectra model (right) of NGC 612 and NGC 3081. Left: The black crosses, red 
filled circles, magenta open circles, and blue crosses represent the data of the FI-XIS, BI-XIS, HXD/PIN, and BAT, respectively, with Icr 
error bars. The spectra of the XIS and PIN are folded with the detector responses in units of counts s — 1 keV -1 , while those of the BAT 
are unfolded spectra in units of photons cm" ~ 2 ks" 1 keV -1 . The best-fit models are plotted by solid curves, and the residuals in units of x 
are shown in the lower panels. Right: The best-fit spectral model in units of EFg (where E is the energy and is the photon spectrum). 
The black, dashed red, dotted blue, dot-dot-dashed cyan, dot-dashed magenta curves correspond to the total, transmitted one, reflection 
component, scattered component, and iron-K emission line, respectively. The purple dashed model below 2 keV in NGC 3081 represents 
the emission from an optically-thin thermal plasma (see text). 
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Table 4 

Best-fit Spectral Parameters with Analytical Models 







NGC 612 


NGC 3081 




Best-fit model 


B 


B ~\~ VcipGC a 




ArGal (i()22 rm ~2\ 

I V jjj \ ^ LI 11 J 


0.0195 


\J . UuOO 


(2) 


M„ fin 22 cm- 2 ) 

1 V J-j^ \ l_^lll ) 


111 ± 5 


98 ± 4 


(J) 




U - &4 -0.07 


n so -i- n 07 
u.oy m u.u < 




r 


1 90 ± 04 


1 co _|_ n no 


(5) 


/scat (%) 


0.55 ± 0.06 


0.73 ±0.07 


(6) 


B cen (keV) 


6 42+ 03 

"•^-0.02 


6.41 ±0.02 


(7) 


E.W. (keV) 


0.28 ± 0.06 


0.35 ±0.05 


(8) 


E.W. refl (keV) 


0.7 ±0.1 


1.0 ±0.1 


(9) 


JVg* (10 22 cm" 2 ) 


(=N H ) 


16t1 


(10) 


R 


0.6 ±0.2 


0.9 ±0.2 


(11) 


F2-10 (ergs cm -2 s _1 ) 


1.6 x 10~ 12 


2.7 x 10" 12 


(12) 


Fio-so (ergs cm -2 s _1 ) 


1.9 x 10" 11 


3.1 x lO" 11 


(13) 


£2-10 (ergss^ 1 ) 


3.0 x 10 43 


3.0 x 10 42 




X 2 /d.o.f. 


100.4/91 


199.3/199 



Note. — (1) The hydrogen c olumn density of Galactic ab- 
sorption by IKalberla et al.l 120051) . (2) The line-of-sight hydro- 
gen column density for the transmitted component. (3) The iron 
abundance relative to the Solar value. (4) The power-law pho- 
ton index. (5) The fraction of the scattered component relative 
to the intrinsic power law. (6) The center energy of the iron- 
K emission line at the rest frame of the source redshift. (7) 
The observed equivalent width of the iron-K line with respect 
to the whole continuum. (8) The observed equivalent width of 
the iron-K line with respect to the reflection component (9) The 
linc-of-sight hydrogen column density for the reflection compo- 
nent. (10) The relative strength of the reflection component to 
the transmitted one, defined as R = S1/2-7T, where f2 is the solid 
angle of the reflector viewed from the nucleus. (11) The observed 
flux in the 2-10 keV band. (12) The observed flux in the 10-50 
keV band. (13) The 2—10 keV intrinsic luminosity corrected for 
the absorption. The errors are 90% confidence limits for a single 
parameter. 

a An additional emission from an optically-thin thermal plasma 
is required, modelled by the vapec code with a temperature of 
kT — 0. 26±0. 02 keV and an emission measure of 1.5X10 63 cm~ 3 . 
The iron abundance is linked to that in the absorber of the trans- 
mitted component (see text). 

scattered component, and iron-K emission line, respec- 
tively. For NGC 3081, the additional soft component is 
also included in purple. 

We find that the fraction of the scattered compo- 
nent with respect to the transmitted one for both of 
NGC 612 and NGC 3081 are fairly small (0.5%-0.8%). 
Also, they are heavily obscured with column densities of 
iVn — 10 24 cm~ 2 , and hence can be regarded as nearly 
"Compton thick" AGNs. Figure [5] shows the correlation 
between the reflection component R and the scattered 
component / sca t, superposed on the same plot presented 
in Paper I. NGC 612 and NGC 3081 are located on a sim- 
ilar position to each other in this plot. In Paper I, the 
authors categorize the observed AGNs into two groups: 
"new type" with R > 0.8 and / sca t < 0.5%, and "clas- 
sical type" with R < 0.8 and / sca t > 0.5%. Since our 
targets are placed just between the two types, it is not 
clear to which "type" these AGNs belong if the intrin- 
sic distribution is indeed distinct. It is also possible that 
the distribution is smooth and they actually represent an 
intermediate class bridging the two types. 

4. TORUS MODEL 

Ikcd a et "ail (|2009T ) performed a set of Monte Carlo sim- 
ulations to calculate a realistic reprocessed emission from 
the torus irradiated by a central source, which is assumed 



to be a cutoff power-law model. In the simulation, they 
assume a 3-dimensional axis-symmetric uniform torus as 
illustrated in Figure [6] It is characterized by the half- 
opening angle f? oa , the inclination angle of torus from an 
observer # mc , the hydrogen column density viewed from 
the equatorial plane N^ q , and the ratio r of r m to r out . 

4.1. Method 

To perfo r m sp ectral fitting on XSPEC with the 
llkeda et al.l (|2009D model, we utilize an atable model 
where the resultant spectra from Monte Carlo calcula- 
tion are stored at grids of the torus parameters in FITS 
files. Here we refer to the results of the continuum repro- 
cessed from the torus ( "torus- reflection component" ) and 
those of iron-K fluorescence line, assuming the geometry 
with r = r in /r out = 0.01. The cutoff energy is fixed at 
£-cut = 360 keV throughout our analysis, which is con- 
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Figure 5. The correlation between the strength of the Compton 
reflection component (R = Q/2n) and the fraction of the scattered 
component (/scat) for our targets (filled circles), superposed on 
the same figure taken from Paper I, where open circle and open 
diamond represent "new type" (R > 0.8 and / SC at ^ 0.5%) and 
"classical type" (R < 0.8 and / sca t 0.5%) AGNs, respectively. 



Obs 




Direct + Disk 

^- Reflection 1 

*- Reflection 2 

" ^ ^ ^ — — — *- Scattered 

Figure 6. Cross -section view of the torus geometry assumed in 
llkeda e t al. (20091). The torus structure is characterized by the 
half-opening angle oa , the inclination angle of torus from an ob- 
server 9i nc , the hydrogen column density viewed from the equato- 
rial plane N^ q , and the ratio of r; n to r ou t- The observed lights 
consist of a transmitted component absorbed by the torus (dash- 
dashed red), a reflection component from the accretion disk ab- 
sorbed by the torus (dash-dashed red), two reflection components 
by the torus; "reflection component 1" absorbed by the torus (dash- 
dot-dotted blue) and "reflection component 2" not absorbed by it 
(dashed orange). We also consider a scattered component by the 
surroinding gas, which is not absorbed by the torus (dash-dash- 
dottcd cyan). 
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sistent with the contraints from the Swift/BAT spectra 
for both targets (> 300 keV). Since the fit is performed 
only in the 0.5-100 keV range where the table model 
is available, the choice of E cut hardly affects our results 
as far as it is higher than ^300 keV. We consider two 
cases for the e lemental abundances; "So lar" abundances 
(as defined bv lAnders fc Grevesselll989() and "Subsolar" 
ones where only those of iron and nickel are set to be 0.5 
times Solar values. Finally, the table files have five free 
parameters; N^ q , oa , 8i nc , the photon index T of the 
incident continuum, and its normalization at 1 keV. 

In the assumed geometry of the torus, the line-of-sight 
hydrogen column density Ah for the transmitted compo- 
nent is related to t hat along the equa torial plane (A^ q ) 
via equation (3) in llkeda et all ([20091 ): 



Ah _ r (cos 9 mc - cos oa ) + sin ( 



N v q (l- r ){rcos0 inc + 



sin ( 



inc — oa ) 

— #oa)} 



(1) 



We introduce torusabs (for the fixed Solar abundances) 
and vtorusabs models (for variable abundances) as lo- 
cal models of XSPEC to represent photoelectric absorp- 
tion of the transmitted component, whose line-of-sight 
column density is related to the torus parameters ac- 
cording to the above equation. In these models we take 
into account Compton scattering processes in addition to 
photoelectric absorption, and hence they can be reliably 
used even for the Compton thick case. We adopt the 
photo electric absorption c ross section by | Verner et al.l 
(fl99l) for consistency with llkeda et al.l (120091). As noted 
by Ikeda et al.l (|2009i ). the cross section bv lVerner et al.l 



(| 19961 ) i s more accurate for energies above 1 keV than 
that by IBaluci nska-Clu irch fc McCammonl (JT992), and 
is nearly equal to the NIST XCOM database^ Since 
the results become physically meaningless for obscured 
AGNs if we obtain inc < oa , we impose the condition 
that #i nc > oa + 1° in the fitting process. 

The fraction of the scattered component to the trans- 
mitted component, / SC at, should be proportional to the 
opening solid angle of the torus if the column density of 
the scattering gas is constant. Thus, we have the con- 
straint that 



cos na = 1 



/sc 



1st 



(1 - COs6»oa i0 ) • 



(2) 



where we normalize / scat = / sca t,o at 8 oa = 9 oa , = 45°. 
Similarly, we also developed the fscat model for XSPEC 
to calculate the normalization of the scattered emission 
as a function of two parameters, 9 oa and / sca t,o- Here 
the normalization parameter / sca t,o reflects the averaged 
column density of the scattering gas and is treated as a 
free parameter. We allow it to vary within 0.1%-5%; note 
that a typical value in Seyfert 2 galaxies is / scat = 3% 
(jGuainazzi et al.1120051 ). 

In addition to the reprocessed emission from the torus, 
we should also expect a reflection component from the 
accretion disk in AGN spectra. Thus, in the transmitted 
compo nent, we include this effect by utilizing the pexrav 
model (Magdziarz & Zdziarski 1995), which is appropri- 
ate to represent the reflection component from semi- 
infinite plane, like that from accretion disks. Here we fix 

12 http: / /www.nist . gov/pml/data/xcom/ index. cf m 



the strength of the disk reflection to be R = Q,/2tt = 1, 
where fi is the solid angle of the accretion disk. The in- 
clination angle of the accretion disk is linked to that of 
the torus. Although the contribution of this component 
is not inclu d ed in the incident photon spectrum in the 
llkeda et al.1 ([20 09) model, the effects on the "torus re- 
flection" spectra are only second-order and are negligible. 
For simplicity, hereafter we refer to the total spectrum in- 
cluding the reflection component from the accretion disk 
as "torus model" . 

To summarize, we can write the torus model of the 
photon spectrum F (E) without the Galactic absorption 
as follows: 

F(E) = 

cxp {-A H (A^ q , 9 oa , inc ) a (E)} I (E) 

/scat.o) I {E) 

+ exp{-A H (n^, 

^oa? ^inc^ O -Rdisk (0inc, E) 

+-Rrefl,l (Ay 13 , oa , iK 

+-Rrefl,2 (A^ q , oa , 9 lr 

+£Fc^Fc (A^ q , oa , 9 U 

+S{E), (3) 

where 1(E) = AE~ T exp (-E / E cut ) is the intrin- 
sic cutoff power-law component, A is the normaliza- 
tion parameter of the intrinsic cutoff power law at 
1 keV, A^ q is the hydrogen column density of the 
torus viewed from the equatorial plane, 9 oa is the half- 
opening angle of the torus, 9 lnc is the inclination an- 
gle of the torus, / sca t.o is the normalization parame- 
ter at 9 oa = 45°, Ah ( N^ 1 , oa , inc J is the absorp- 
tion column density for the transmitted component, 
a (E) is the cross section of photoelectric absorption, 
fscat iPoai fsc&t, o) is the scattered fraction, i? disk (0 inc , E) 
is the Compton reflection component from the accre- 
tion disk, i? r on,i (A^ q , oa , 0inc) is the torus-reflection 

component 1 (see Figure EJ), i? ro fi,2 (a^, oa , inc ) is 

the torus-reflection component 2 (see Figure O, £Fe is 
the normalization parameter for the iron-K emission 

line, Lpc ^A^ q , oa , 0j nc ^ is the iron-K emission line, and 

S (E) represents additional soft components (for the case 
of NGC 3081 where the apec model is used). We allow 
the relative normalization of the iron-K line £fo with re- 
spect to the torus-reflection components to float between 
0.5 and 2, in order to cover the iron abundance range over 
the fixed values (0.5 or 1.0) and to model the effects of 
time variability between the transmitted component and 
the averaged reprocessed emission. 

In short, Equation ([3]) is expressed as absorbed trans- 
mission^] + scattcrin£3 + absorbed accretion disk re- 



13 In XSPEC nomenclature, torusabs*zhighect*zpowerlw 

14 In XSPEC nomenclature, fscat *zhighect*zpowerlw 
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Figure 7. The observed spectra (left) and the best-fit spectral model (right) with the torus model. Left: same as Figure [4] (left) . Right: 
the best-fit spectral model in units of EFg (where E is the energy and Fg is the photon spectrum); total (thick black), transmitted 
component (thick dashed red), reflection component from the accretion disk (thin dot-dashed green), torus reflection component 1 (thin 
blue), torus reflection component 2 (thin dotted orange), scattered component (thin dot-dot-dashed cyan), iron-K emission line (thin 
dot-dashed magenta). The purple dashed curve below 2 keV in NGC 3081 represents the emission from an optically-thin thermal plasma. 



flectiorP^I + absorbed torus reflectiorPI + vmabsorbed 
torus reflectiorFl + iron lin<H There are eight free pa- 



find that the column density in the equatorial plane, 



rameters: £> inc , oa , N^ q , / scat ,o, ^ccn, £fc, A, and T. 
In the spectral fit, we employ the MINUIT MIGRAD 
method ("migrad") as the fitting algorithm, which is 
found to be more stable than the standard Levennerg- 
Marquardt method ("leven") in our case utilizing the 
numerical models. Both of the Suzaku and BAT spectra 
(but below 100 keV) are used throughout this section. 

4.2. Application 
4.2.1. NGC 612 

Figure [7] plots the best- fit torus model, whose param- 
eters are summarized in Table [5] The torus model with 
"Subsolar" abundances is adopted for this target, based 
on the fitting result with the analytical models (Sec- 
tion 3). We find that 6> oa ~ 60° - 70° and 6> inc > 76°. As 
shown in Figure [7j the contribution from the reflection 
component 2 (unabsorbed one) is very small, suggest- 
ing that we are seeing the target from an edge-on angle. 
In such case, the inclination angle 9i nc can be poorly 
determined above a certain threshold, because the ob- 
served line-of-sight hydrogen column density is rather 
insensitive to #j nc for a given N^ q value. Actually, we 

15 In XSPEC nomenclature, torusabs*toruspexrav 

16 In XSPEC nomenclature, atable{refll_torus.fits} 

17 In XSPEC nomenclature, atable{refl2_torus.flts} 

18 In XSPEC nomenclature, con- 
stant *atable{refl_fe_torus. fits} 
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24.1 



cm , is close to that along the line of 



sight (ATh) as estimated from the analytical model fit. 
We also find that / sca t,o — 0.14%, which indicates that 
the amount of scattering gas around the nucleus is re- 
markably small. 

4.2.2. NGC 3081 

The best-fit parameters with the torus model for 
NGC 3081 are summarized in Table [5] and the model is 
plotted in Figure[7l For this target, we adopt the "Solar" 
abundance tables based on the analytical model fit. We 
obtain 6> oa ~ 15°, 8 inc ~ 19°, iV"| q ~ 10 240 cm~ 2 , and 
/ SC at,o - 4.6%. By contrast to NGC 612, / scat ,o is rather 
large, and there is a significant contribution from the 
unabsorbed reflection component, indicating that the in- 
clination angle must be close to the torus opening angle. 
The results are consistent with the picture proposed for 
the "new type" AGN SWIF T J0601.9-8636 (ESO 005- 
G004) bv lUeda et "all (|2007t ). that the nucleus is deeply 
buried in the geometrically thick torus and is observed 
from a rather face-on angle. 

5. DISCUSSION AND CONCLUSION 

With Suzaku follow-up, we have obtained the best- 
quality broad band spectra covering the 0.5-60 keV band 
of two Swift/BAT AGNs, NGC 612 and NGC 3081. 
First, we found a range in the iron abundance; NGC 612 
has about 0.5 times Solar abundance of iron (where "So- 
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Table 5 

Best-fit Spectral Parameters with Torus Model 



NGC 612 NGC 3081 



(1) 


Table model 


Subsolar 


Solar + apec a 


(2) 


jyGal ( 1Q 22 cm ~2) 

7v| q (10 22 cm" 2 ) 


0.0195 


0.0388 


(3) 


U3±J° 




(4) 


# oa b (degrees) 


70 (> 58) 


15 ±2 


(5) 


6> inc (degrees) 


87 (> 76) 


19±X 


(fi) 


r 


1.9 ±0.1 


2.0 ±0.1 


(7) 


/scat.O (%) 


0.14 (< 0.19) 


4.6 (> 3.1) 


(8) 


Been (keV) 


6 4S+ - 11 


6 42+ ' 03 


(9) 






0.57 (< 0.69) 




X 2 /d.o.f. 


84.3/83 


209.0/200 


Note. 


— (1) The tabic model used in the fit. 


"Solar" means 



the table with Solar abundances, while "Subsolar" means the 
table with 0.5-time iron and nickel abundances with respect to 
Solar o nes. (2) The hydrogen column density of Galactic absorp- 
tion by Kalbcrla ct al. (2005). (3) The hydrogen column density 
of the torus viewed from the equatorial direction. (4) The half 
opening angle of the torus. (5) The inclination angle of the torus. 
(6) The power-law photon index. (7) The fraction of the scat- 
tered component relative to the intrinsic power law when the half 
opening angle of the torus is 45°. (8) The center energy of the 
iron-K emission line at the rest frame of the source redshift. (9) 
The relative strength of the iron-K emission line to that predicted 
by the torus model. The errors are 90% confidence limits for a 
single parameter. 

a An additional emission from an optically-thin thermal plasma 
with Solar abundances is required, modelled by the apec code 
with a temperature of kT — 0.25 ± 0.02 keV and an emission 
measure of 1.5 X 10 63 cm -3 (see text). 

k The range of the 6 oa , is limited to < 70° in the torus model. 

lax" corresponds to Fe/H = 4.68 x 10~ 5 ), which is sig- 
nificantly smaller than that of NGC 3081. Applying 
the analytical models, we find that these objects are 
nearly Compton thick AGNs with Vh — 10 24 cm -2 and 
the fraction of the scattered component with respect 
to the transmitted component is small, / sca t < 0.8%, 
suggesting that these be l ong to "hidden" population ac- 
cording to I Winter et al.l (|2009f ). Plotting the results in 
the /scat versus reflection strength (R) plane, we find 
these two targets are located just between those occu- 
pied by "new type" (geometrically thick tori) and "clas- 
sical type" AGNs defined in Paper I, implying that they 
would be an intermediate class bridging the two types. 
We need a larger sample to reveal the true distribution 
of the whole AGN population in this plane. In this 
context, simultaneous broad band observations of more 
"new type" candidates are important to examine their 
reflection strengths, such as those with small scatter- 
ing fractions identifi ed from the XMM-Newton catalog 
(jNoguchi et al.ll2009h . 

To further investigate the details of the torus geome- 
try of the two AGNs, we apply numerical spectral mod- 
els based on Monte Carlo simulation where a simple 3- 
dimensional geometry of the to rus i s assumed, fo l lowin g 
the work bv llkeda et ail (j2009ft and lAwaki et al.l (|2009t ). 
We also consider the Compton reflection component from 
the accretion disk. To our knowledge, this is the first 
time all effects both from the torus and disk are self- 
consistcntly considered in spectral analysis of obscured 
AGNs. It is remarkable that we are able to reproduce 
the observed spectra quite well with this torus model, 
which has only 3 free geometrical parameters; the open- 
ing angle, inclination, and equatorial column density. 



The column density along the equator plane is found 
to be Vh ~ 10 24 cm -2 for both sources, which is also 
si milar to that found from the Seyfert 2 galaxy Mrk 3 
by llkeda et al.l <|2009f ) . The relative absence of higher 
column densities, though very limited in number, may be 
consistent with the fact that even hard X-rays above 10 
keV have a bias against detecting heavily Compton thick 
AGNs with Vh ~ 10 25 cm" 2 , unless the sample i s limited 
to the very local universe (Mali zia et all 12009ft . Thus, 
a majority of Swift/BAT AGNs do not have extremely 
Compton thick tori defined at the equator plane unless 
observed from a face-on angle. Future sensitive hard X- 
ray surveys may start to pick up such populations, whose 
number density and cosmological evolution are still open 
questions. 

The analysis with the torus model suggest that the 
torus geometry of the two targets may be different in 
spite of the very similar results obtained from the ana- 
lytical models. Our results confirm that the fundamental 
assumption of the unified model where the opening an- 
gles are all the same is too simple. For NGC 612, we find 
that the opening angle is relatively large (~ 60° — 70°) 
and the object is observed from an edgc-on angle, con- 
sistent with a picture of "classical type" Seyfert 2 galax- 
ies. Similar torus parameters are obtained for Mrk 3 by 
llkeda et all (|2009f) . By contrast, the torus opening angle 
of NGC 3081 is much smaller (~ 15°), and we observe 
it from a face-on angle. This implies that N GC 3081 is 
closer to a "new type" AGN discovered by lUeda et al.l 
( 2007) surrounded by a geometrically thick torus. This 
picture for NGC 3081 is consistent with the time vari- 
ability of the column density, because we are seeing the 
thinnest part of the torus t hat is expected to be highly 
patchy (Risal iti et al.l 12002). We note, however, that 
the best-fit torus parameters we obtain from the present 
analysis should not be taken at their face values, which 
could depend on the initial assumpt ion of the torus ge - 
ometry. For instance, as discussed in llkeda et al.l ([2009) , 
if we assume a lower value for r(= r; n /r ou t) than 0.01, 
we would obtain a slightly larger half-opening angle 9 oa 
for the same N^ q and 9 lnc to account for the increased 
contribution of the unabsorbed reflection component. 

Since the observed scattering fraction is similar be- 
tween the two targets, this difference in the torus open- 
ing angle indicates that the amount of scattering gas 
around the nucleus is much smaller in NGC 621 than 
in NGC 3081, as represented in the obtained / sca t o 
value, / sca t o < 0.2% for NGC 612 and / scat > 3% for 
NGC 3081.' The small amount of the gas in NGC 612 may 
be consistent with its classification as a "weak emission 
line" radio galaxy, where the jets expel the surrounding 
gas. By contrast, the detection of the optically-thin com- 
ponents in NGC 3081 could represent the abundance of 
the ambient gas around the nucleus. 

An important implication from the present study is 
that the classification of different types of tori (e.g., ge- 
ometrically thin or thick) based solely on the scattered 
fraction may be difficult in some cases. Our work has 
demonstrated the power of the application of numerical 
torus models based on Monte Carlo simulation to best 
extract the physical view of the nucleus beyond the sim- 
ple phenomenological spectral analysis, although caution 
must be paid because we have considered only the sim- 
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plest geometry by assuming a uniform density. Combina- 
tions of the high quality broad band X-ray spectra with 
more realistic numerical simulations will be a key ap- 
proach for further understanding of the nature of AGNs. 
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